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ABSTRACT 

An  extension  of  a boundary  integral  method  without  using  Green’s  function  for  solving 
electromagnetic  boundary-value  problems  in  layered  media  is  presented.  The 
generalization  of  the  previous  applications  of  that  method  to  the  case  under 
consideration  is  achieved  by  using  fractional  cylindrical  functions  as  the  testing 
functions,  which  account  for  spatial  inhomogeneity  of  the  ambient  medium.  To 
illustrate  advantages  of  such  approach,  numerical  analyses  are  presented  for  2D 
scattering  problem  involving  penetrable  cylinder  of  elliptical  cross-section  shape. 

INTRODUCTION 

Employment  of  fractional  cylindrical  functions  for  analysis  of  electromagnetic  wave 
scattering  in  the  presence  of  bodies  with  coordinate  surfaces  was  offered  in  [1,2].  But  in 
our  opinion,  such  functions  can  be  used  to  solve  many  others  problems  of  modern 
electromagnetic  theory.  Therefore,  the  attempt  of  expanding  implementation  area  of 
such  functions  was  undertaken.  Two-dimensional  scattering  problem  for  homogeneous 
cylinders  of  arbitrary  cross-section  shape  embedded  in  a plane  layer  is  considered. 

PROBLEM  STATEMENT  AND  BASIC  EQUATIONS 

Let  suppose  that  as  it  shown  in  Fig.l  a permeable  cylinder  of  arbitrary  cross  section  S p 

is  situated  in  one  layer  of  the  three-layer  structure.  Introduce  a coordinate  system  XYZ 
and  suppose,  that  the  impress  sources  f of  monochromatic  (~exp  (-ico  !))  wave  and 
layers  boundaries  belongs  to  surrounding  domain  Se . Cylinder  formative  is  parallel  to 
axis  OX.  The  medium  inside  inclusion  is  described  by  material  parameter  e (F)  and 


Fig.  1.  Geometry  of  the  problem. 
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wavenumber  kp  (F) , and  outside  of  inclusion  - zc(r)  and  k(r)  correspondingly. 
Fields  amplitudes  outside  and  inside  of  inclusion  satisfies  to  such  equations: 


e(?)Vx~Vi+*2(F) 
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They  also  satisfy  continuity  conditions  for  functions  ue , 1/e,  duJdN,  u , 1/e, 

dup/dN  on  structural  boundaries  in  Se,S  and  ue , and  satisfies  radiation  condition 

when  r ->  +oo . N - normal  to  L - contour  of  the  scatterer,  directed  to  Se . In  the  case  of 

excitation  by  vertically  polarized  wave  (h  = {//.,0,0},i?  = {o,^,#.}),  it  is  suitable  to 

choose  an  x-  components  of  magnetic  field  vectors  as  a functions  ue,  « . As  a result  of 

manipulations  described  in  [3],  one  can  obtain  following  formulae  for  the  scattered 
field: 


w = f_L 
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where  G{r,r')-  Green's  function  of  a regular  problem.  We’ll  represent  an  unknown 
functions  ue  and  duJdN  as  an  expansions  in  terms  of  functions  and  d^c/dN  : 

«,(F)  = ^a0^(F);  due(r)/dN  = Mr)/dN  . (6) 

a a 

Unlike  previous  papers,  we'll  choose  the  function  £,a(F)  in  a following  way: 

/„,  U)  cos(mcp)  a = (m,e),  m = 0,1,2,.. 

£,a(F)  = < ; cp  = arg(F) , (8) 

f„,  (r)sin(mcp)  a = (m,  o),  m = 1,2,.. 


where  fm(r)  = Jm+\(r)/ Jm(r) , Jm(r ) - Bessel  functions  of  order  m.  We’ll  compute  the 
wave  functions  using  the  fractional  cylindrical  ones,  as  it  was  shown  in  [3]. 

A possibility  of  fractional  cylindrical  functions  application  to  the  considered  problem 
justified  by  presence  of  cylindrical  inclusion  with  sufficiently  smooth  contour  of  cross- 
section.  Unlike  the  results  obtained  in  [2],  the  discussed  problem  is  complicated  by 
presence  of  dielectric  layer  and  cylinder  with  sufficiently  arbitrary  cross-section. 


DESCRIPTION  OF  NUMERICAL  ALGORITHM  AND  ANALYSIS  OF 
OBTAINED  RESULTS 

In  this  section  some  numerical  results  based  both  on  the  analytical  techniques  developed 
in  previous  papers  and  a new  algorithm  for  evaluating  scattering  diagrams  by  means  of 
fractional  cylindrical  functions  are  presented.  The  goal  of  numerical  experiments  was  to 
determine  effectiveness  of  the  proposed  approach.  Computations  performed  for 
dielectric  elliptic  cylinders,  immersed  in  homogeneous  dielectric  layer.  They  show  that 
for  cylinders  with  transversal  dimensions  less  than  wavelength  application  of  fractional 
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cylindrical  functions  does  not  give  noticeable  advantage.  This  fact  can  be  easily 
explained,  if  one  take  into  account,  that  the  basic  advantages  of  fractional  functions 
related  to  their  smooth  behavior  attached  to  large  index  and  argument  values  [2].  Also 
it's  necessary  to  mention,  that  the  most  effective  this  procedure  became  in  a case  of 
circular  cylinder. 

In  Fig.  2 scattering  diagrams,  calculated  by  means  of  Bessel  and  fractional  cylindrical 
functions  are  shown.  Dielectric  permeabilities  of  cylinder,  layer,  upper  and  lower  half- 
spaces are  equal  to  ep  =1.2,  es  =3.0.  cc  =1.0.  st  =2.0  accordingly.  Cylinder’s  radius 

is  2.0 -X,  layer’s  thickness  - h-6T.  an  inclusion  embedding  depth  - Zp=3  C. 

Incident  wave  length  - X = 1 .5  m.  One  can  see.  that  when  Bessel  functions  are  used, 
scattering  diagrams  became  steady  only  when  M>30  (M  maximum  order  of  the 
scattering  matrix).  On  the  other  hand,  for  fractional  cylindrical  functions  analogous 
results  can  be  obtained  when  M=20.  and  additional  increasing  of  number  of  terms  in 
decomposition  does  not  leads  to  any  changes. 


Observation  angle  <t> 

a) 


b) 


Fig. 2. Scattered  field  diagram  obtained  for  Bessel  functions  (a)  and  for  fractional  cylindrical  functions  (b). 


CONCLUSION 

Obtained  results  show,  that  the  fractional  cylindrical  functions  within  the  framework  of 
the  null  field  method  forms  a successive  base  for  solving  a variety  of  scattering 
problems.  Such  problems  in  particular  are  essential  elements  in  microwave  and  optical 
device  design,  nondestructive  testing,  remote  sensing  and  thin  film  physics.  The  same 
technique  can  be  applied  to  another  diffraction  problems  (for  example  3D  problems) 
with  only  minor  modifications. 
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